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We employ a combination of the ab initio band structure methods and dynamical mean-field 
theory to determine the electronic structure and phase stability of paramagnetic FeO at high pres¬ 
sure and temperature. Our results reveal a high-spin to low-spin transition within the B1 crystal 
structure of FeO upon compression of the lattice volume above 73 GPa. The spin-state transition 
is accompanied by an orbital-selective Mott metal-insulator transition (MIT). The lattice volume 
is found to collapse by about 8.5 % at the MIT, implying a complex interplay between electronic 
and lattice degrees of freedom. Our results for the electronic structure and lattice properties are in 
overall good agreement with experimental data. 

PACS numbers: 71.10.-w, 71.27.-|-a, 71.30.-|-h 


I. INTRODUCTION 

Iron monoxide, wiistite (FeO), is a basic oxide compo¬ 
nent of the Earth’s interior. Its electronic state and phase 
stability is of fundamental importance for understanding 
the properties and evolution of the Earth’s lower mantle 
and outer core.— Nevertheless, in spite of long-term inten¬ 
sive research, the phase diagram of FeO at high pressures 
and temperatures, as well as several key properties, such 
as electronic structure and magnetic properties, are still 
poorly understood. 

FeO has a relatively complex pressure-temperature 
phase diagram with at least five allotrops.— Under am¬ 
bient conditions, it is a paramagnetic Mott insulator 
with a rock-salt B1 crystal structure. Upon compres¬ 
sion above ~ 16 GPa, Bl-type FeO undergoes a struc¬ 
tural transition to the rhombohedral i?3 phase (rBl),— 
which further transforms to the NiAs (B8) structure 
above 90 GPa— Under ambient pressure, FeO is anti¬ 
ferromagnetic below the Neel temperature Tjv 198 K.— 
The high-pressure properties of FeO have attracted much 
recent interest both from a theoretical and experimen¬ 
tal point of views— Shock-wave compression and elec¬ 
trical conductivity experiments reported a possible ex¬ 
istence of a high-pressure metallic phase of FeO above 
~ 70 GPa— ^ On the basis of high-pressure Mossbauer 
spectroscopy measurements?^ the metallic state was as¬ 
signed to a high-spin (HS) to low-spin (LS) transition, 
which has been proposed by first-principles calculations 
to occur above ^ 100-200 GPa— In contrast to that, 
x-ray emission spectroscopy indicates that the Fe high- 
spin state is preserved at least up to ^ 140 GPa (at 
room temperature))^ while collapsing to the LS state 
upon further heating—^ On the basis of these measure¬ 
ments, the insulator-to-metal HS-LS transition has long 
been considered to be due to a structural transformation 
from the B1 to B8 lattice. Nevertheless, recent exper¬ 
iments have shown that the Bl-type FeO undergoes a 
high-temperature insulator-to-metal transition at about 
70 GPa, retaining the same lattice structureJ^ The ob¬ 
served metallic phase was proposed to relate to a spin- 


state crossover. Moreover, it has been shown that the 
Bl-type structure remains stable at high pressure and 
temperature, being the stable phase along the geotherm 
through the Earth’s mantle and outer core— *22 

These recent experiments have lead us to reinvestigate 
the electronic structure and local magnetic state of Ee 
in EeO at high pressure and temperature by employ¬ 
ing a combination of the ab initio electronic structure 
methods and dynamical mean-field theory of strongly 
correlated electrons (LDA-I-DMET)— (LDA stands 
here for the local density approximation). Applications 
of LDA-I-DMET have shown to capture all generic as¬ 
pects of a Mott metal-insulator transition, such as a co¬ 
herent quasiparticle behavior, formation of the lower- 
and upper-Hubbard bands, and strong renormalization 
of the effective electron mass, providing a good quan¬ 
titative description of the electronic and lattice prop¬ 
erties. In particular, we employ an important advance 
of the LDA-t-DMFT approach which is able to deter¬ 
mine the electronic structure and phase stability of cor¬ 
related materials We use this advanced theory to 
investigate the electronic structure and phase stability 
of paramagnetic FeO at high pressure and temperature, 
which remained unexplored up to now. We find that 
magnetic collapse occurs in the B1 crystal structure of 
paramagnetic FeO upon compression to ^ 73 GPa, in 
agreement with experiment. The HS-LS transition is 
intimately linked with an orbital-selective Mott metal- 
insulator transition, which is accompanied with a col¬ 
lapse of the lattice volume by about 8.5 % at the MIT. 
Our results for the electronic structure and lattice prop¬ 
erties are in overall good agreement with experimental 
data. 


II. METHOD 

In this work, we investigate the electronic structure 
and phase stability of paramagnetic FeO under pres¬ 
sure using the GGA-f-DMFT computational approach 
(GGA: generalized gradient approximation). To this end. 
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we calculate the total energy and local moments of the 
B1 cubic crystal structure of FeO as a function of lat¬ 
tice volume. Below we denote the compressed phase 
by the relative volume w.r.t. the calculated equilib¬ 
rium lattice volume as v = I^/Vq. We employ a fully 
charge self-consistent GGA-I-DMFT approach^ imple¬ 
mented with plane-wave pseudopotentials^. For the par¬ 
tially filled Fe 3d and O 2p orbitals we construct a ba¬ 
sis set of atomic-centered symmetry-constrained Wannier 
functions^i^. To solve the realistic many-body problem, 
we employ the continuous-time hybridization-expansion 
quantum Monte-Carlo algorithm^. The calculations 
are performed in the paramagnetic state at temperature 
T = 1160 K. We use the average Coulomb interaction 
U = 7 eV and Hund’s exchange J = 0.89 eV for the Fe 
3d shell, in accordance with previous estimatesid^ii^iii. 
The Coulomb interaction is treated in the density-density 
approximation. The spin-orbit coupling is neglected in 
these calculations. The U and J values are assumed 
to remain constant upon variation of the lattice vol¬ 
ume. In addition, we check how our results depend on 
a possible reduction of the U value under pressure. For 
this purpose, we perform calculations with a substan¬ 
tially smaller value, [7 = 5 eV. We employ the fully- 
localized double-counting correction, evaluated from the 
self-consistently determined local occupancies, to ac¬ 
count for the electronic interactions already described 
by GCA. The spectral functions were computed using 
the maximum entropy method. The angle resolved spec¬ 
tra were evaluated from analytic continuation of the self¬ 
energy using Fade approximants. 
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FIG. 1: Top: Total energy and local moment of para¬ 

magnetic FeO calculated by GGA-I-DMFT as a function of 
lattice volume. The lattice collapse associated with the HS- 
LS state transition is depicted by a vertical black dashed line. 
Bottom: Fe 3d and partial t 2 g/eg occupations as a function 
of volume. 


III. RESULTS AND DISCUSSION 

We first compute the electronic structure of FeO within 
non-spin-polarized GCA using the plane-wave pseudopo¬ 
tential approaches Overall, our results agree well with 
previous band-structure calculations. In particular, the 
calculated equilibrium lattice constant is found to be a = 
7.74 a.u. The calculated bulk modulus is B = 232 GPa. 
We note however that GGA calculations give a metal¬ 
lic solution, in qualitative disagreement with experiment. 
The calculated equilibrium lattice constant is remarkably 
smaller, by 6-7 %, than the experimental one. We also 
note a strong overestimation of the bulk modulus, which 
is more than 25 % off the experimetal value. Glearly, 
standard band-structure techniques cannot explain the 
properties of paramagnetic FeO, since they do not treat 
electronic correlations adequately. 

To resolve this obstacle, we now compute the electronic 
structure and phase stability of FeO using the fully charge 
self-consistent GGA-l-DMFT method. In Fig. [T] (top) we 
show the e voluti on of the total energy and local magnetic 
moment \/(wfy of paramagnetic FeO as a function of lat¬ 
tice volume. We fit the calculated total energy using the 
third-order Birch-Murnaghan equation of states^ sepa¬ 
rately for the low- and high-volume regions. Overall, our 


results for the electronic structure and lattice properties 
of FeO, which now include the effect of electronic corre¬ 
lations, agree well with experimental data. In particular, 
our calculations at ambient pressure give a Mott insulat¬ 
ing solution with an energy gap of ^ 0.8 eV. The energy 
gap lies between the top of the valence band originating 
from the mixed Fe 3d and O 2p states and empty Fe 4s 
states. Our result for the Mott d-d energy gap is about 
2 eV, in good agreement with optical and photoemission 
experiments!^ We hnd the equilibrium lattice constant 
a = 8.36 a.u., which is less than 1-2 % off the experi¬ 
mental value. The calculated bulk modul us is B = 140 
GPa, the local magnetic moment ^ 3.7 ps- Fe 

t 2 g and eg orbital occupancies are 0.68 and 0.55, respec¬ 
tively. These findings clearly indicate that at ambient 
pressure Fe^+ ion is in a high-spin state (^ = 2). In fact, 
in a cubic crystal field, the Fe^'*' ions (i.e., Fe 3d^ conhg- 
uration with four electrons in the t 2 g and two in the Cg 
orbitals) have a local moment of 4 pB- 

Furthermore, the total-energy and local-moment calcu¬ 
lation results exhibit a remarkable anomaly upon com¬ 
pression of the lattice volume down to v ^ 0.72. In¬ 
deed, the local moment is seen to retain its high-spin 
value down to about v 0.8. Upon further compres¬ 
sion, a broad high-spin (HS) to low-spin (LS) crossover 
takes place, with a collapse of the local moment to a 
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FIG. 2: Evolution of the spectral function of paramagnetic 
FeO as a function of lattice volume. Fe t 2 g/eg and O p orbital 
contributions are shown. The metal-insulator transition asso¬ 
ciated with the HS-LS transformation takes place at i/ ~ 0.72, 
at pressure ~ 73 GPa. 


LS state with magnetic moment ~ 1.2 ps at pressure 
above 160 GPa, i.e., v < 0.6. We display our results for 
the evolution of the Fe tgg and Cg orbital occupations in 
Fig. [T] (bottom). Upon compression, we observe a sub¬ 
stantial redistribution of charge between the t 2 g and Cg 
orbitals within the Fe 3d shell. Fe t 2 g orbital occupa¬ 
tions are found to gradually increase with pressure, re¬ 
sulting in (almost) completely occupied state (tgg occu¬ 
pation is about 0.95). On the other hand, the Cg orbitals 
are strongly depopulated (their occupation is 0.25) and 
the Fe 3d total occupancy remains essentially unchanged 
with pressure. We therefore interpret this spin crossover 
as a HS-LS transition. 

The spin-state transition is accompanied by a struc¬ 
tural transformation. The structural change takes place 
upon a compression of the lattice volume to v 0.72, 
resulting in a collapse of the lattice volume by ^ 8.5 %. 
This value should be considered as an upper-bound es¬ 
timate because we neglect multiple intermediate-phase 
transitions when fit the total-energy result to the third- 
order Birch-Murnaghan equation of states. Our estimate 
for the transition pressure is p = 73 GPa. In addition, we 
note that the bulk modulus in the HS phase is somewhat 
smaller than that in the LS phase (162 GPa), imply¬ 
ing an enhancement of the compressibility at the phase 
transition. Overall, the electronic structure, the equi¬ 
librium lattice constant, and the structural phase sta¬ 
bility of paramagnetic FeO obtained by the fully charge 
self-consistent GGA-l-DMFT approach are in remarkably 
good agreement with the experimental data. Our find¬ 
ings clearly indicate the crucial importance of electronic 
correlations to explain the electronic structure and lattice 
properties of paramagnetic FeO. 

Next we address the spectral properties of paramag¬ 
netic FeO. In Fig. [2] we present the evolution of the 
spectral function of FeO as a function of lattice vol¬ 
ume. We note that at ambient pressure paramagnetic 


FeO is a Mott insulator with a relatively large d-d gap 
~ 2 eV. Interestingly, the energy gap lies between the 
Fe t 2 g states, while the t 2 g-eg gap is somewhat larger, 
~ 2.5 eV. The top of the valence band has a mixed Fe 
3d and O 2p character, with a resonant peak in the filled 
t 2 g band located at about -0.9 eV. The latter can be as¬ 
cribed to the formation of a Zhang-Rice bound state4^ 
We find the d-d energy gap being gradually decreased 
upon compression, resulting in an orbital-selective Mott 
metal-insulator transition (MIT). In fact, the onset of 
the local moment collapse at ^ ~ 0.75 is associated with 
closing of the gap for the t 2 g orbitals, while the Cg states 
remain to be gaped down to ly ~ 0.7. We note that 
the HS-LS transition takes place upon compression of 
the lattice volume by v ~ 0.72, at critical pressure of 
p ~ 73 GPa, in agreement with experiment. Our calcu¬ 
lations clearly show that the HS-LS state transition in 
the HI structure of paramagnetic FeO is associated with 
an orbital-selective Mott-Hubbard MIT. Moreover, the 
MIT is accompanied by a collapse of the lattice volume 
by about 8.5 %. 

In addition, we check how the electronic structure of 
FeO depend on a (possible) reduction of the Coulomb in¬ 
teraction U with pressure. For this purpose, we perform 
calculations with a substantially smaller value U = 5 eV 
and the same Hund’s coupling value J = 0.89 eV. 
Our new results are in overall qualitative agreement 
with those presented above. We notice that the HS- 
LS crossover is accompanied by an orbital-selective MIT 
which is found to occur at somewhat smaller compres¬ 
sion V 0.77, at critical pressure 55 GPa. We also 
note a substantial reduction of the Mott d-d band gap 
to about 1.1 eV, whereas the charge-transfer gap be¬ 
tween the Fe 3d-0 2p and empty Fe 4s states remains 
unchanged, about 0.8 eV. These hndings clearly indicate 
that our results for the spin-state crossover and associ¬ 
ated with it the orbital-selective MIT are robust. The 
HS to LS transition takes place even at a substantially 
reduced (possibly due to the applied pressure) value of 
U. 

We also calculated the momentum-resolved spectral 
function of paramagnetic FeO. In Fig. [3] we present our 
results obtained for the HI structure of paramagnetic 
FeO across the MIT. Our calculations at ambient pres¬ 
sure show a Mott insulator with an energy gap value of 
^ 0.8 eV [see Fig. [3] (bottom)]. The energy gap lies be¬ 
tween the occupied states with a mixed Fe 3d-0 p char¬ 
acter and empty Fe 4s state (the latter is clearly seen 
as a broad parabolic-like band at the F-point just above 
the Fermi energy), in agreement with photoemission and 
optical experiments4i In fact, the optical spectroscopy 
measurements observe a weak absorption between 0.5 
and 2.0 eV, assigned to the Fe 3d-0 2p to Fe 4s tran¬ 
sitions. The strong absorption edge associated with the 
d-d transitions is found to appear in optical spectroscopy 
at about 2.4 eV. Our result for the d-d gap at ambient 
pressure is about 2 and 2.5 eV for the t 2 g-t 2 g and t 2 g- 
Cg transitions, respectively. The O 2p states are about 
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FIG. 3: The k-resolved spectral function of paramagnetic FeO 
computed within GGA+DMFT for different volumes. Top: 
Results for the LS phase with lattice volume v ~ 0.6. Bottom: 
The equilibrium HS phase. 


-4 eV below the Fermi level, but still have a substan¬ 
tial contribution near the Fermi level. We note however 
that the Fe tig states are seen to have a predominant 
contribution around the Fermi level. We therefore in¬ 
terpret paramagnetic FeO as a Mott-Hubbard insulator, 
in agreement with previous studies. Furthermore, our 
calculations show an entire reconstruction of the elec¬ 
tronic structure of paramagnetic FeO in the LS phase 
[see Fig. [3] (top)]. Namely, we obtain a strongly corre¬ 
lated metal with almost completely occupied tig band, 
which is located in the vicinity of the Fermi level, and 
a broad disperse spectral weight originating from the Fe 
Cg band crossing the Fermi level. Moreover, the Fe tig 
states exhibit a resonant state just below the Fermi level 
which can be ascribed to the Zhang-Rice bound state. 
The O 2p states are shifted to about -8 eV below the 
Fermi level, and now have a moderate contribution near 
the Fermi level. 

We also calculate the local (dynamical) spin-spin cor¬ 
relation function x(r) = {rhz(T)rhz{0)) of paramagnetic 
FeO for different volumes, where r is the imaginary time. 
In Fig. 0] we display our results for the corresponding 
intra-orbital tig and eg contributions. We note that the 
tig contributions are seen to be almost independent of 
T, even in the LS state, implying that the tig moment 
remains localized across the HS-LS transition. On the 



FIG. 4: Local spin-spin correlation function xi’’') calculated 
by GGA+DMFT for paramagnetic FeO as a function of vol¬ 
ume. The intra-orbital tig and Cg contributions are shown. 



V/V„ 

FIG. 5: Quasiparticle weight Z and spectral weight at the 
Fermi level iV(e_F) = —^G{t = 13/2) calculated for the Fe 
tig and Cg orbitals across the MIT in paramagnetic FeO. The 
structural change associated with the HS-LS transition is in¬ 
dicated by a vertical black dashed line. The orbital-selective 
MIT phase is marked by a red filled rectangle. 


other hand, the Cg states clearly exhibit a crossover from 
localized to itinerant magnetic behavior under pressure, 
suggesting that the HS-LS metal-insulator transition in 
paramagnetic FeO is of the orbital-selective type. In ad¬ 
dition, we notice that at high pressures the screened local 
moment defined as oc dTx{T) differs from the cor¬ 
responding instantaneous (m^) value (mostly because of 
the contribution originating from the Cg band), which 
suggests an itinerant-moment behavior of the LS state. 

Finally, we calculate the quasiparticle weight employ¬ 
ing a polynomial fit of the imaginary part of the self¬ 
energy E(ia;„) at the lowest Matsubara frequencies 
It is evaluated as Z = [1 — dImTj{iuj )/from the 
slope of the polynomial fit at a; = 0. In Fig. [S] we present 
our results for the partial tig and Cg contributions com¬ 
puted for different volumes. In the LS phase, the tig and 
eg quasiparticle weights are finite and found to decrease 
monotonously upon expansion of the lattice. Upon fur¬ 
ther expansion of the lattice, the tig and Cg Z factors are 
seen to collapse to zero at different volumes. This result 
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again suggests that the HS-LS transition is accompanied 
by an orbital-selective MIT. The electronic effective mass 
diverges at the MIT, in accordance with a Brinkman-Rice 
picture of the MIT4i We note that this divergence coin¬ 
cides with the drop of the spectral weight for the t 2 g and 
eg orbitals at the Fermi level shown in Fig. [S] 

IV. CONCLUSION 

We determined the electronic structure and phase sta¬ 
bility of the B1 crystal structure of paramagnetic FeO at 
high pressure and temperature. Our results clearly estab¬ 
lish a high-spin to low-spin transition which is found to 
appear upon compression of the lattice by U/Vq 0.72, 
at pressure 73 GPa, in agreement with experiment. The 
spin-state transition is intimately linked with an orbital- 
selective Mott-Hubbard metal-insulator transition. The 


MIT is accompanied by a collapse of the lattice volume 
by ^ 8.5 %, implying a complex interplay between elec¬ 
tronic and lattice degrees of freedom. Upon compres¬ 
sion, the Fe Cg states exhibit a crossover from a localized 
to itinerant-moment behavior, while the t 2 g moment re¬ 
mains localized across the HS-LS transition. Our results 
are consistent with the picture of an electronically driven 
spin-state transition and volume collapse in the BI crys¬ 
tal structure of paramagnetic FeO. 
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